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Abstract 

We consider a SIRS model with time dependent transmission rate. 

We assume time dependent vaccination which confers the same immu¬ 
nity as natural infection. We study two types of vaccination strate¬ 
gies: i) optimal vaccination, in the sense that it minimizes the effort 
of vaccination in the set of vaccination strategies for which, for any 
sufficiently small perturbation of the disease free state, the number of 
infectious individuals is monotonically decreasing; ii) Nash-equilibria 
strategies where all individuals simultaneously minimize the joint risk 
of vaccination versus the risk of the disease. The former case corre¬ 
sponds to an optimal solution for mandatory vaccinations, while the 
second corresponds to the equilibrium to be expected if vaccination 
is fully voluntary. We are able to show the existence of both optimal 
and Nash strategies in a general setting. In general, these strategies 
will not be functions but Radon measures. For specific forms of the 
transmission rate, we provide explicit formulas for the optimal and 
the Nash vaccination strategies. 
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1 Introduction 


Vaccination is the best response available in the control of most infectious 
diseases. A huge effort is put on the development of new and better vaccines. 
When humans are directly involved, the role of direct experimentation is nat¬ 
urally limited and therefore mathematical models have been used to evaluate 
the effect of control measures, such as vaccination, to assist in policy deci¬ 
sions. One central result of classical mathematical models for the spread of 
infectious diseases is that persistence of an infectious disease within a pop¬ 
ulation requires the density of susceptible individuals to exceed a strictly 
positive critical value such that, on average, each primary case of infection 
generates more than one secondary case. It is therefore not necessary to 
vaccinate everyone within a community to eliminate infection. This phe¬ 
nomenon is known as herd immunity and is one of the key epidemiological 
questions in defining a vaccination strategy. 

In this work, we consider a SIRS model with periodic transmission. The 
model consists of a non-autonomous system of ordinary differential equations 
in which we introduce periodic vaccination of adults. For simplicity, we 
considered that vaccination confers the same protection as natural infection. 
We study the consequences of two extreme types of vaccination strategies: 
mandatory vaccination, where the population is vaccinated at a predefined 
rate; and voluntary vaccination, where individuals can choose freely to be 
vaccinated or not, according to their risk perception. 

Mathematical models have been widely used to help health authorities 
in the definition of vaccination strategies for very different contexts. Typi¬ 
cally, the objective is to define an optimal vaccination strategy by minimizing 
combinations of the vaccination effort/cost and of the effective reproduction 
number IZq (i.e., the number of secondary infections generated by a primary 


case) (Muller and Hadeler, 

1996 

Castillo-Chavez and Feng, 

1998 

Laguzet 

and Turinici, 2015a 

). This 

can give rise to particularly interesting problems 


when we consider non-homogeneous models for which vaccination strategy 


depends on age (I 

duller and Hadeler, 1996 C 

lastillo-Chavez and Feng, 

1998 

Tartof et al. 

20K 

1), risk-groups (jScott et al. 

2015; 

Long and Owens 

2011) 

or when it is time-dependent ( 

Onyango and Muller 

, 2014 d’Onofrio 

2002 

Browne et al. 

2015; Houy 20H 

). The current work is concerned with time- 


dependent epidemic models with vaccination when both the transmission 
rate and the vaccination are assumed to be periodic. Note that periodic 
vaccinations are used by public health services; e.g., influenza vaccine is only 
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available in a specific season of the year. 

Whenever the goal is long term disease elimination, optimal vaccination 
will consist on reducing IZo below one, in the sense that it implies the at¬ 
tractiveness and the asymptotic stability of the disease-free state. Here, we 
choose to work with an alternative definition of optimal vaccination, where 
the goal is, not only to eliminate disease but also to prevent outbreaks. 
Hence, we define a class of preventive vaccination profiles such that, for any 
sufficiently small perturbation of the disease free state, the number of in¬ 
fectious individuals is monotonically decreasing. We construct the optimal 
vaccination strategy as the limit of preventive strategies for which vaccination 
effort is minimized. We start by revisiting basic concepts in mathematical 
epidemiology to recall that for constant transmission, the condition TZq < 1, 
and subsequent stability of the disease-free state, is equivalent to the condi¬ 
tion that infectious population / is monotonically decreasing in time, when 
the initial number of susceptible individuals is below the number of suscep¬ 
tible individuals in the disease free state. However, as we move towards 
more general situations this equivalence may not hold. Note that the for¬ 
mer condition refers to the long term behaviour of the system, while the 
latter considers also the short time behaviour which, in principle, is more re¬ 
strictive (Hastings, 2010, 2004). In particular, the former condition restricts 
the average transmission over a period and the latter is defined pointwise in 
time. Our approach is particularly suitable for diseases with high mortality 
or morbidity rates, for which it is imperative to prevent outbreaks. Despite 
that, so far our model and examples will not consider disease related death. 

On the opposite end of vaccination policies is voluntary vaccination, which 
is increasingly common in industrialized countries. Even when vaccines are 
offered by the public health system without costs, vaccination is, at least in 
part, voluntary. Opposition to vaccines can be philosophical, religious and 
depend on social contacts and information available. It puts important chal¬ 
lenges to disease control by decreasing vaccine uptake. The case which is best 
known is the measles, since the unproven hypotheses that measles-mumps- 
rubella (MMR) vaccine was linked to autism led to a decrease in vaccination 


followed by measles epidemics in UK (Fitzpatrick, 2004 Jansen et ah, 2003) 


Voluntary vaccination can also give rise to free-rider phenomenon, where in¬ 
dividuals or families choose not to be vaccinated, or to not have their children 
vaccinated, taking advantage of herd immunity created in the population by 
others, avoiding the possible negative effects of vaccination. In this work, 
we consider a population of rational individuals that compares the risk of 
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the vaccination (more precisely, its perception of the risk of the vaccination) 
and the risk of the disease and make options that minimizes the joint risk. 
Despite the fact that some countries are implementing hires for parents that 
prefer not to vaccinate their childrerQ, we do not introduce in the model a 
risk of non-vaccination, other than the one associated to the disease. 

In this work, we model human behaviour using game theory. In a seminal 


paper by Bauch and Earn (2004), it has been shown that voluntary vaccina¬ 
tion cannot lead to disease eradication. The authors coupled a SIR model 
for disease spread in a partially vaccinated population with a theoretical 
game framework describing a population of rational individuals. Many sub¬ 
sequent developments were made in order to include the human behaviour in 


epidemiological models (cf. Chen 

(2006) 

; d’Onofrio et al. (i 

>007); 

Manfredi 

et al. (2009) 

Coelho and Codego 

(2009); 

Mbah et al. (2012 

; Manfredi and 

D’Onofrio (2 

013); f 

dor in et al. 

(21 

113); Bhattacharyya et al. 

(2015 

; Laguzet 

and Turinici 

(20151 

:>)). See also (1 

bank et al., 2010; 

Wang et ah, 2 

015) for a 

review, and 

Funk et ah, 2015 

for further discussion on the subject. 


In this paper, we generalize the framework of Bauch and Earn (2004) to 


the SIRS model with periodic transmission function. From the modelling 
point of view, we consider that all choices in the population influence the 
dynamics, and the resulting dynamics also has effect in the rational behaviour 
of the population. Due to the richness of the non-autonomous system that 
describes our model, several technical problems arise. For instance, the risk 
of disease no longer depends simply on the constant steady state as before. 
Considering a rational individual, we assume that he/she is going to choose 
to be vaccinated only when the risk of disease times the probability of being 
infected is higher than the risk of the vaccine, as perceived by the taker. 
As we analyse only stationary states of our periodic system, the risk to 
be minimized is the joint risk of vaccination and disease during one season. 
Hence, we define the set of herd immunity provider vaccination strategies, for 
which the rational strategy for a given focal (rational) individual is not to be 
vaccinated, taking advantage of the herd immunity provided by the choices of 
the rest of the population. Moreover, we define a Nash vaccination strategy 
as the strategy that minimizes the joint risk for every individual taking into 
account the strategy of all other individuals, i.e., the natural strategy to be 


1 That’s the case of Poland and Australia. See http://www.thenews.p1/l/9/Artykul/204007, 
Parents-fined-for-not-vaccinating-children and http://naturalsociety.com/australia- 
enforces-15k-penalty-for-parents-who-dont-vaccinate/, respectively. 
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expected in a population of rational individuals with full knowledge of all 
epidemiological data. 

Existence of optimal and Nash vaccination strategies are proved in this 
work in a very general setting; however, these strategies may not be functions 
but Radon measures, even for transmission rates given by a real function. 
This is a consequence of the fact that the set of continuous functions in a 
given compact interval is not closed under any reasonable metric. Many 
results used in the existence proofs presented in the appendices require com¬ 
pactness and after introducing a convenient topology in the set of continuous 
functions, we are naturally led to the introduction, in this framework, of 
Radon measures. For more information on the topic of Radon measures we 
refer to (Schwartz 1973| Athreya and Lahiri 2006). 

The paper is organized as follows. In Section [2j we introduce the math¬ 
ematical model and derive some preliminary results. Section [3] is dedicated 
to the vaccination strategies. First, we give rigorous definitions of preventive 
vaccination strategies, and of vaccination effort and we define the optimal 
strategy as one strategy that can be arbitrarily approximated by a preven¬ 
tive strategy and such that its associated effort is never superior to the effort 
of any given preventive strategy. In the context of voluntary vaccination, we 
define the set of herd immunity provider strategies and the concept of Nash 
strategy, in which all individuals minimizes the joint risk of vaccination and 


disease. In the end of the Subsection 3.2, we state the main theorem, which 
guarantees the existence of an optimal and a Nash vaccination strategies in 
the set of Radon measures. Explicit formulas for the optimal and Nash strate¬ 


gies are provided in Subsect ion |3.3[ for specific forms of the transmission rate. 
In Section [4j we present some examples such as the constant transmission 
case, the sinusoidal case and also a critical case to illustrate the results from 
previous sections. We finish with two appendices, the first one guaranteeing 
the existence of periodic solutions in the model and the second proving the 
existence of optimal and Nash strategies. 


2 The Model 


Consider a SIRS model. Let S(t), I(t), R{t ) be the fraction of susceptible, 
infectious and recovered individuals at time t > 0. We assume non negative 
normalized initial conditions, i.e, S'(O),/(0), A(0) > 0, S'(0)+/(0) + R(0) = 1. 

/3 „ T T 7 T~> T~> a O O P 


We also assume the transitions S + I 


21, I 


A, R —>■ S, S 


R. 
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Figure 1: Schematic diagram of the SIRS model with death and vaccinations. 


Constants p (mortality/birth rate), a (temporary immunity) and 7 (recovery 
rate) are strictly positive. The disease is assumed to be non-fatal, i.e., the 
death rate p > 0 does not depend on the disease class. By normalization, 
we also consider the birth rate as / 1 . These are common assumptions of 
compartmental epidemiological models. 

We consider functions (3,p : M + —> M + , representing the transmission and 
vaccination rates at time t, respectively. More precise assumptions on these 
functions will be introduced latter on. 

From now on, we call SIRS model to the following system of differential 
equations: 


/j + aR — f3(t)IS — p(t)S — pS 


p + a — al — /3(t)IS — p(t)S — (p + a)S 

(la) 

/3(t)IS — 7 / — pi 

(lb) 

7 / + p(t)S — pR — aR 

(lc) 


A schematic representation of the SIRS model with vaccination is represented 
inFigure[l| Due to the normalization S(t)+I(t)+R(t) = S'(0)+/(0)+R(0) = 
1, the equation for R is always redundant and will be ignored from now on. 
We define A 2 := {(x, y ) e M 2 | x > 0, y > 0, x + y < 1}. 


We begin by analysing the solutions and stability of system ([Tj) . 

Lemma 1 . Let us consider that functions (3 and p are continuous functions 
with commensurable periods, i.e., there exists T > 0 such that p(t + T) = p(t ) 
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and j3(t + T) = (3{t ) for all t. Equivalently, we assume that /3,p G (7(10, T]) 
with (p(0),/3(0)) = (p(T),/3(T )). 

Therefore, there exists only one periodic solution of system ([Tj) in the 
subspace {/(t) = 0, \/t}, given by (5 0 (f), 0) G A 2 . We call this solution the 
disease-free solution. This solution attracts all initial conditions of the form 
(5(0),0). We define I 0 (t) = 0. 

Depending on the choices of the parameters 7 , p and the functions p , (3, 
we may have one of two possibilities: 

1. The disease-free solution is globally stable in A 2 ; 

2. There are other periodic solutions (with period multiple ofT), called 
endemic solutions (Sift), lift)) G A 2 , with lift ) > 0, for all t G M+ and 
i G N. In this case, there is 7 > 0 such that for any initial condition 
(5(0),/(0)) with /(O) > 0, we have liminf^oo Iff) > 77 . In this case, 
we say that the solution of the SIRS model is persistent. 


Furthermore, for any initial condition, the solution (S(t),I(t)) depends con¬ 
tinuously on (3 and p; namely, if p n —* p and (3 n —>■ (3 weakly as measures 
and both sequences are uniformly integrable, then Si\p n ,/3 m ] —>■ 5j[p, (3\ and 
Ii\p n ,/3 m ] —> Ii(p, (3] uniformly in [0, T\, when n,m —> 00 , i G N. 


Proof. For the first part, see Rebelo et al. (2012); see Appendix [A] for further 
details. The second part follows from Heunis (1984, theorem 2.1). □ 


Remark 1. The assumptions on (3,p in the previous lemma are extremely 
restrictive and used only for the first part of the result (existence of dis¬ 
ease free solution mid periodicity of the endemic solution). If we relax our 
assumptions to require only that (3 is of bounded variation and p is a mea¬ 
surable function, then existence of solutions (not necessarily periodic) and 
convergence of solutions (as in the second part of Lemma [7]) is guaranteed 
by (Heunis. 198f). This will be explored in the examples. Note that if (3 n ,p n 
are continuous then, (3 := lim/3 n and p := limp n necessarily satisfy these 
more relaxed assumptions. 


For constant transmission and vaccination, we establish the following re¬ 
sult, that is going to motivate our definition of optimal vaccination. 

Lemma 2. Consider that j3(t ) = (3 0 > 0 and p(t ) = Po > 0. The only sta¬ 
tionary disease free solution of system 0 is given by 5 0 = 5 0 [p] = 5 0 [p](t) = 
f +a , < 1. Furthermore, the three conditions below are equivalent: 

po+U+a — ’ ^ 
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^ 0qSq < 7 
7+M — 


2. The disease free solution (S 0 , 0) is globally asymptotically stable. 

3. I' < 0 for all I > 0 and all S < S 0 . 


Proof. 


1^3 


Note that I' = 1(7 + p) (777 — lj, and therefore, assuming 

if anrl nnlir if Q ^ 7 


I > 0, I' < 0 if and only if S' < Sp 


1 
3 

EN2T 1 


ho 

S < S 0 < 777 and then I' < 0. 

I' < 0 if and only if S' < and therefore S 0 < Sj(r- 


ho 


[Ml We follow ideas from ( ]Cruz[ |2009[ ) ; see also ( |Capasso| |1993[ ) for 


other examples of use of Lyapunov functions in Mathematical Epidemiology. 
Let us define 

Z(s,/) = i(s -7 + /)%<>/, a = 2 7±i!I±7±jy 

We differentiate V with respect to t and obtain 
V'{S, I) = ^ = (s - So + /) (S' + /') + aV 


— (S — S 0 + IJ ((So — S)(p 0 + p + ci) — (7 + ci + p)/) + a/' 

= — (s — S 0 j^ (po + p + cx) — / 2 (y + ol + p) 

— (S — So)/(2(p + a) + 7 + po) + «(So-^S — (7 + p))/. 

= — (S — SqJ (po + p + 0 ) — J 2 (7 + ci + p) 

'7 + P 


a/3oI 


So 


Sr 


Let G = {(S,/) G [0,1] x (0,1]|S + / < 1}. Note that V' is a continuous 
function in G and, by Condition [lj it is negative, ft is immediate to verify 
that V is a Lyapunov function in G. Let G be the closure of G. As {(S, I) G 
G\V'(S,I) = 0 } is the s ingleton with the equilibrium point { (So, 0 )} we 
(2009, Corollary 1.2 in Chapter X) that (So, 0) is globally 


conclude from 


asymptotically stable. 


Hale 






















HMD After system linearisation around the disease free solution 
(Sq, 0), we find the Jacobian matrix 


-Po - h 
0 


a 


—a — (3o§o 


If > 1, the equilibrium would not be stable, which leads to a contradic¬ 
tion. 

0 

From the above lemma, we recover the effective reproductive number for 
the constant parameter case, 1Z 0 := 7 Z 0 [p\ : = . Condition 1Z 0 < 1 

guarantees at the same time that all epidemics will be eventually extinct and 
that /(f) decreases monotonically in time, from 7(0) > 0. 

However, in the time dependent case (in particular in the periodic case), 
these two phenomena are not equivalent. In general, even for linear systems, 
it is possible that before being attracted to an asymptotic equilibrium, the 
trajectory of (5(f),/(f)) drifts away from this equilibrium (HastingsJ 2010, 
2004). 

For the periodic case, we can compute the effective reproduction number, 


following (Thieme, 2000) (see also (Wang and Zhao, 2008)), as 


Tin : — 


7 + h 


(PS 0 ) = 


T (7 + h) Jo 


/5(f)5 0 (f)df . 


( 2 ) 


Note that, for the periodic case, condition 1Z 0 < 1 still guarantees asymp¬ 


totic stability of the disease free case (Wang and Zhao, 2008), but does not 


necessarily prevent the existence of outbreaks; see, e.g., (Zhao, 2008). 


In this work, we will look for conditions that generalize, for time-dependent 
parameters, Condition [3] in Lemma[2j i.e., that guarantees that the number of 
infectious is monotonically decreasing for small perturbations of the disease 
free solution. From the modelling point of view, no particular definition can 
be considered better than the other; in fact, for certain particular diseases 
(e.g., polio, tuberculosis) vaccination policy aims to eradicate/eliminate the 
disease in the long run, while for other diseases, governments act to prevent 
the existence of large outbreaks (e.g., influenza, cholera) (WHO, 2015). Our 
approach describes better this second setting. 
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From now on, we assume that, for a given vaccination strategy p(t), sys¬ 
tem I(t )) is in its stationary (periodic) state, and we will consider two 

different cases: 

(Cl) The disease free state (<So[p] (t), 0); 

(C2) A certain endemic state (Si[p](£), Ii[p](t)). (There is no uniqueness 
for the endemic state; for the sake of simplicity, we will consider from 
now on only one endemic solution. There is no essential change if we 
consider more than one.) 

Both solutions are assumed to be periodic, possibly with period multiple of 
T; however, without loss of generality, we will consider the period given by 
T. Note that for a different set of parameters more complicated behaviour 
(possible chaotic) can be found, cf. (Kuznetsov and Piccardi, 1994). 


3 Vaccination strategies 

In this section we will consider two types of vaccination: mandatory and 
voluntary vaccination. For each one, we will define one special case: for the 
former, an optimal vaccination is defined as one vaccination strategy that is 
able to prevent outbreaks while having the minimum number of vaccinations 
possible and for the latter, a Nash vaccination strategy is defined as a strategy 
in which all individuals in population minimize the joint risk of both disease 
and vaccine. 


(Kuznetsov and Piccardi, 1994 


3.1 Optimal vaccination 


For the optimal vaccination, we choose to work with a generalization of 
Condition [3] in Lemma [2] More specifically, we say that a certain vaccination 
strategy p is a preventive strategy when the fraction of individuals in the 
class / decreases monotonically in time for any small enough perturbation of 
the disease free state. We then construct the optimal vaccination strategy 
as the limit of the preventive strategies for which the vaccination effort is 
minimized. In (Onyango and Muller, 2014), optimality for time dependent 


vaccination profiles is defined based on the effective reproductive number. 
Note that in our model, only susceptibles are vaccinated, which implies a full 
knowledge of the current status of an individual. 
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Definition 1. We define the vaccination effort associated to a given strategy 
p, as the average number of vaccinations in one period, i.e. E[p] := (pS[p\) := 
Jq p(t)S[p](t)dt, where (•) denotes the average in one period, and S is the 
relevant solution, given by (T|7J) or <®. defined above. 

We denote a cumulative distribution function, associated with p, by 
P(i) = f 0 dp, or in a more relaxed notation dP(i) = p(t)dt. To simplify 
the notation, we will use indistinctly dP and pdt, whenever there is no risk 
of confusion. Therefore, we now write E[p] = y / 0 S'(r)dP(r). For technical 
reasons, we need to consider bounds in the set of vaccination profiles. More 
precisely: 


Definition 2. We say that a certain vaccination function p is admissible if 
its cumulative distribution is such that 

P([0,T])= f dP<( A + a )jW ( 3 ) 

Jo 7 + h 

where f3 = sup t6 r 0 T i /3(t). Furthermore, we use RM + to denote the set of non¬ 
negative Radon measures P in [0, T\ such that P([0, T]) <(p + a an d 

C + ([0, T]) to denote the set of continuous functions in [0, T), withpiff) = p(T ) 
such that fy p(t)dt < (p + a)T^J^. We also consider the natural immersion 
C + ([0,T])cRM + . 


Now we show that the definition of vaccination effort can be extend con¬ 
tinuously for the case of Radon measures. 


Lemma 3. Let p n e C+([0, T]) be such that p n —>■ p G RM +; 
topology, cf. (Koralov and Sinai, 2001), and lefE\p] := limE[p„] 
is independent of the choice of the sequence p n . 


in the weak 
Then, E[p] 


Proof. Let p n ,q n G C'+QO, T]) such that p n ,q n —> p G RM + . Let P„ and Q n 
be the associated cumulative distribution functions, respectively. Note that 


E \p n \ ~ E[g n ] 



S[q n ] dQn) 

Qn) + [ T (S[p n \-S[q n ])dQ n 


From the fact that S\p n ] is bounded and P n — Q n —>■ 0, we conclude that 
the first integral converges to 0. For the second integral, the convergence 
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to 0 follows from the continuity of p H- S[p] in the appropriate topology. 
See (Heunis, 1984, theorem 2 . 1 ) for further details. O 


Definition 3. Let f3 6 C([0,T]) ; /3(0) = (3(T) be given. For a given vacci¬ 
nation strategy p, let (5o(t),0) = (/So[p](t), 0) be the disease-free solution of 
System 0 - We say that p is a preventive strategy if (3(t)So[p\(t) < 7 + p 
for all t. We call Xp = Xp[fi\ the set of admissible strategies that are preven¬ 
tive, i.e, x P = {p e C+QO,! 1 ]), with p( 0) — p(T), and /3(t)S 0 [p](t) <7 + p, 
for all t e [0, T]}. 


Now, we analyse the preventive strategies. First, we explicitly charac¬ 
terize the disease free state and then we show the existence of at least one 
preventive strategy. Afterwards, we define the concept of optimal strategy. 
Here, we reproduce the result from (Thieme 2003, Theorem 3.7). 


Lemma 4. Let S 0 [p] (t) be the time dependent periodic number of susceptibles 
in the unique disease free state of system (JT|) . Then 


So[p](0) 


(p + a) f Q T e -/ s T (p+^+“)(Hdr ds 
l — e - If (p+/i+a)(r)dr 


Before looking for optimal strategies, we prove that the set of preventive 
strategies is not empty. 


Lemma 5. For any choice of parameters, there exists at least one preventive 
strategy, i.e., xM 7 ^ 0 f or all fd G C([0, T]), f3( 0) = P(T). 


(J-- 

\7+M 


constant, corre- 


Proof. Given (3, assume p(t) = p 0 > (p + a) 
sponding to a preventive strategy in the case of the maximum transmission 
rate. Note that Sq\jj] (t) = v ^+“ +Cl is the only stationary solution of the Equa¬ 
tion (la) with I{t) = 0 and is, additionally, the solution obtained from the 


initial condition given by Lemma [4] From the definition of p^, we conclude 
that (3(t)S 0 \p](t) < 7 + p for all t, and then {p(t) = p 0 } G Xp- □ 


Finally, we construct the optimal vaccination strategy as one strategy that 
can be arbitrarily approximated by a preventive strategy and such that its 
associated effort is never superior to the effort of any preventive strategy. 
More rigorously, we define an optimal vaccination strategy by 
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Definition 4. Let (3 G C([0, T\), /3(0) = /3(T) be given. We say that a given 
strategy p opt = Pop t [/5] is optimal if the following conditions are simultaneously 
satisfied: 


1. there is at least one sequence x P [/3] 9 p n —> p opt (in measure). 

2. for any p G x p [/3], E [p\ > E[p opt \. 

Whenever p opt is a function (as discussed in Subsection 3.3), Definition [4] 
means simply that there is a sequence of preventive strategies that converge 
to p 0 pt and that p opt minimizes the vaccination effort in the closure of the set 
X P - 


3.2 Rational vaccination 

In this subsection, we study a population of rational individuals and how 
their decisions influence the disease dynamics. 

For each focal individual the probability of getting the disease is assumed 
to depend on the disease incidence for each time step. The following lemma 
shows how this probability can be computed from the model. 

Lemma 6. The probability that a susceptible non-vaccinated mdividual at 
time t gets the disease between times t and t + At, for At sufficiently small, 
is given by (3(t)I(t)At + o (At). 

Proof. All susceptible non-vaccinated individuals are in the category S. From 
time t to time t + At, (3(t)I(t)S(t)At individuals will be infected, pS(t)At 
will die and the remainder S(t) — [pS(t) + (3(t) I(t)S(t)] At will be in the class 
S at time t + At. Therefore, the probability to be infected from times t to 
t + At, given that he/she did not die, is given by 

__ = 8(t)I(t)At + o (At) 

S(t)-[yiS(t)+f3(t)I(t)S(t)]At + /3(t)I(t)S(t)At ,uu 1 + 

where we used (1 — pA(f )) -1 = 1 + o (At). 

See Figure [ 2 j for a schematic representation of this reasoning. □ 

For the voluntary vaccination, we consider that a rational individual will 
(not) vaccinate him/herself if the risk of the disease times the probability to 
get the disease, given the overall strategy of the population, is larger than 
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Figure 2: Transitions of non-vaccinated individuals from state S at time t 
out by death (down arrow) and to states S and I at time t + At. Note that 
indications in the arrows are for the total number of individuals leaving state 
S during interval At. 

(respectively, small than) the risk of the vaccine. If both risks are the same, 
any strategy is equally advantageous. A fully informed rational individual 
will access, in the beginning of the season, the probability to get the disease, 
using all available epidemiological data, and decides his/her personal strategy 
as the strategy that minimizes the joint risk, i.e., the risk of the disease times 
the probability to get it (conditional to no vaccination), plus the risk of the 
vaccine (conditional to vaccination), during the next season. 

We start by defining the set of immunity provider strategies, i.e., the 
set of strategies for which a focal rational individual will decide to be not 
vaccinated. 

Definition 5. Let (3 G (7([0, T]), /3(0) = (3(T ) be given. For a given vaccina¬ 
tion strategy p, assume the existence of a persistent endemic solution (Si, Ji). 
Let r d > 0 and r v > 0 be the risks of the disease and of the vaccination, re¬ 
spectively. We define r := —. We say that p is a herd immunity provider 
strategy if /3(t)Ii\p](t) < r for all t. We call Xh = Xh[(3\ the set of all herd 
immunity provider strategies, i.e. 

Xh = {p £ C'+QO, T]) with p( 0) =p(T), and f3(t)Ii\p\(t) < r,for all t e [0, T]} . 

If there is no endemic solution, we define Xh — {p€C + ([0, T]) with p(0) —p(T)j 

Note that, from the definition, it is clear that any preventive strategy is 
also herd immunity provider, i.e., Xp C Xh- 
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Finally, we will define the Nash-equilibrium strategy as the strategy that 
minimizes the joint risk for every individual, given the strategy of all other 
individuals. 

Definition 6 . Let (3 G (7([0, T]), (3(0) = (3(T ) be given. Let us consider a 
population with strategy p G C+([0, T]), and a focal individual that uses vac¬ 
cination strategy p* G RM + . Let P and P* be the cumulative distributions, 
associated to p and p *, respectively. Assume that the focal individual is sus¬ 
ceptible at time t = 0, and therefore the probability to be susceptible at a later 
time t is given by e~Jo The joint (disease arid vaccination) risk during 

one season (i.e, the probability that something bad — disease or reaction to 
the vaccine — happens in one season, times the associated risks) is given by 



o Jo 



Given a strategy p, a rational individual will choose a strategy p * such that 
for every strategy p' G RM + 


p\p*,p\ < p\p',p\ ■ 


We say that pNash £ RM + is a Nash strategy if for any sequence p n G 
C + ([0, T]), such that p n —> pNash and for every strategy p' G RM + , 


limsup (p[p N ashj Pn\ ~ P[P , Pn]) < 0. 


If PNash is a function, the above definition simplifies to the assertion that 
p[pNash,PNash] < p[p',PNash] for every strategy p' G RM + . 

We finish the subsection stating the existence theorem for both optimal 
and Nash-equilibrium strategies. In general terms, for (3 G C([0,T]), with 
0(0) = (3(T), we prove that there is at least one optimal vaccination strategy 
and at least one Nash vaccination strategy. These strategies may not be 
functions, but measures. This implies that, after rewriting System (|TJ) in 
the form X' = T(f, A"), the function T : M + x A 2 —> M 2 is a Charatheodory 
function (i.e., measurable in the first variable and continuous in the second) 
and therefore there is a (weakest) topology which guarantees existence of 
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solutions of the differential equations and gives continuous dependence for 
each initial data point. See (Aliprantis and Border 1999 Heunis 1984) for 


further details. The proof of the existence theorem below will be postponed 
to Appendix |B| 


Theorem 1. Assume (3 G C([0,T]), with (3( 0) = (3(T). Then, there is at 
least one optimal vaccination strategy p opt [(3] and at least one Nash vaccina¬ 
tion strategy p Nas h [/5] - 


3.3 Vaccination strategies for regular transmission func¬ 
tions 


Despite the fact that we cannot guarantee a priori existence of optimal and 
Nash strategies as functions, we will provide precise conditions for which p opt 
and/or pNash are functions. In particular, we derive explicit formulas for the 
optimal and Nash strategies for sufficiently regular transmission functions 
j3, with some extra technical conditions. In the end, we discuss vaccination 
strategies when (3 is discontinuous (in particular of bounded variation). 

We start by finding an explicit formula for p opt in some special cases. 


Theorem 2. Consider system ([!]). Assume that 
f3'(t) > -{/i + a)/3(t) 

Then 


Po P t(t) = (n + a) 
is an optimal strategy. 


m 

7 + A 


T —1 

1 

(4) 

7 + h / 

- 1 ) + m 

) m 

(5) 


Proof. First note that p opt (t) > 0 if and only if (3(t ) satisfies Equation Q 
for all t > 0 . 

We divide the proof in several steps: 

1 st step: We start by using Lemma[4]to show that /3(0)S'[p op t](0) =7 + pi. 
Indeed, let p = p opt and therefore 


SoboptKO) 


iogf£? ds 

1 _ e -^Jo T h(r)dT 
rT d — /3(r)dr i 

Jo Ts e 7+A1 dg = 7 + /^ 

l_ e -f$f/o T /3(Ddr ,0(0) 
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2 nd step: Now, we show that for any t > 0, /3(t)S 0 [p op t}(t) =7 + p. Using 
Equation (la) with S(t) = S\jj] (t), I(t ) = /[p](0) = 0 and p = p opt , we find 


(mm' 

= F(t)S(t) + l3(t)S'(t) 


= /3'(t)S(t) + /3(t)(p + a) 


(3(t)(p + a) 


- (3(t)(p + a)S(t ) 


= P(t)(p + a) 


msm 

7 + h / ' 


m 

7 + h 


1 )+m 


S(t ) 


We conclude that f3(t)S(t ) = 7 + p is the unique solution of the last equation 
with the initial condition found in the first step. 

3 rd step: Let P, be the cumulative distribution associated to pi, i = 1,2. 
We will prove now that if [J dP] > dP 2 , s £ [0, T) and / 0 T dPi > J 0 T dP 2 , 
then So[pi](£) < So[p 2 ](t). For simplicity, we write Si = So [pi] , i — 1,2. From 
Lemma |4j it is clear that S'i(O) < ^(O). Furthermore, 

(Si - S 2 )' + Pl (t)(S 1 - S 2 ) + (/i + a)(S' 1 - S 2 ) = — (pi -p- 2 )S 2 < 0 . 


After rewriting the last equation, we find that 
d r 


d t 


J 0 ( P 1 +M+a)(r) dr( Si _ g ^ = _ e / 0 (pi+/.+«)(r)dr^ i _ < 0> 


and conclude that Si(t) < S 2 (t) for all t. 

4 th step: We now show that E[p opt ] = (p + a)(l — (S 0 [p O pt]))- Indeed, let 
p n be a sequence such that f dP n > f s dP opt for any s > 0, where P„ and 
P opt are the cumulative distributions associated to p n and p opt , respectively. 
Assume, furthermore, that p n -£ p opt as measure. We use S n = So[p n ]: 
therefore S n (t) < So[p op t](^) = and then p n £ y p . Furthermore, 0 = 

(S’ n ) = p + a - E \p n \ - (p + a)(S n ) and E [p n ] = (p + a)(l - (S 0 [p n ])). We 
take n — >• 00 and use the continuity of So in p n . This finishes this step. 

5 th step: Finally, we prove that for any p £ y p , E[p] > E[p opt ]. From 
S 0 [p](t) < ^ = S 0 [p op t](t), we conclude that {S 0 [p}) < (S 0 [p O pt]) and there¬ 
fore E[p] > (/i + a)(l - (S 0 [p opt ])) = E[p opt ]. □ 

Now, we show that for a non-constant seasonal epidemics, it is always 
better to consider the natural fluctuations, also at the level of the vaccination 


17 









campaign. This result goes along with Agur et al. (1993). See also the 
discussion in Onyango and Muller (2014). 


Corollary 1. Let f3 be a non-constant periodic function, a7id assume that 
the optimal strategy p op t[l3] is given by Equation 0- Then, E[p opt [/3]] < 
E[Po P t [</?>]] ■ 

Proof. We use the classical harmonic/arithmetic mean inequality, i.e., < 

( (3 ), with equality if and only if f3 is constant. Hence, 


E [Popt \P}] = (/i + a) (1 - (So [popt [/?]])) = (fi + a)(l- (^^)) 

<(v + <x)(i- = (p + «) (! - s ° b°pt K^)]]) 

= E [popt [<£>]], 


with equality if and only if [3 is constant. □ 

Corollary 2. In the optimal vaccination case, the total number of vacci¬ 
nations in a single season do not exceed the number of newborns plus the 
number of individuals that lost their immunity during the previous season, 
i.e, assuming the worst case scenario, where everyone got the disease or was 
vaccinated, i.e., E[p op t] < p + a. 

Proof. Consider po the vaccination strategy defined in Lemma [5] Then 

E[p op t] < E[po] = k f Po(t)S 0 [p 0 }(t)dt = P() ^ < p +a . 

j- JO Po + p + 01 


□ 

The preceding corollary shows that it is possible to prevent outbreaks 
without vaccinating the entire population, which means that it is still possi¬ 
ble to attain the herd immunity effect in this more restricted framework of 
preventive strategies. 

We finish this section by stating an explicit formula for pNash in some 
special cases. 
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Theorem 3. Assume r < inf^^r] P{t), 


m 

m 


< 7 + H < 


m 

m 


+ P(t) 


and 


d 
d f 


e (r+a- 7 )t 


d 

df 


0 ( 7 +v)t \" 

7 w) 


< e (r+n+a)t 


(p, + a - 



Then, the strategy given by 


( 6 ) 

(7) 


PNash(f) 


m 


(7 + T)P{t)-/3'(t) [ 

— [r + p + a] 


is a Nash-equilibrium strategy. 


s P'(t) (f'(t)) 2 P"(t) 


ar 


nt) nt) v 2 (t) m 


( 8 ) 


Proof. We show that p Nas h > 0, Si(t),h(t) e [0,1], Vf and /3(f)/i[p N ash](f) = 
r for all t > 0. Initially, let us define 


Sl (t) = -?M= e -(7 


£(*) /3 2 (t) 


5 (7+//)t 


dt y /3(f) 


Note that Equation ([6]) guarantees that Si(f) e [0,1] for all t. Also, the 
assumption on r implies that /(f) G [0,1] for all t. Furthermore, 

.. p + a - S[ (f) - al At) . . 

PNash(t) = ---— - (r + p + a) . 


With this definition, note that 


-/3(t)SJi = S[ +p N ash5'i - (p + a)(l - Si) + ah = -rSi 

and therefore /3(f)/i[pNash](f) = r > 0. Finally, we use Equation (J7|) to prove 
that pNash > 0. From Definition [bj we conclude that p[p*,PNash] does not 
depend on p* and therefore PNash is a Nash-equilibrium. O 

Note that the definition of Si was motivated by plugging /(f) = r//3(t ) 
into Equation ( |lb| ); in particular, this choice of /(f) cancels the first term of 
p on Definition [6j 
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3.4 Comparisons 

Now, we compare the vaccination effort associated with the two extreme 
vaccination strategies. 

Proposition 1. Assume that /J^ash ft Xp■ Then, E[pNash] < E[p opt ]. 

Proof. We use that 0 = (/3IiSi) — (7 + /r)(/i) and 0 — (n + a)(l — (Si)) — 
E[p N ash] - a(h) ~ {PhSi) to conclude that E[p Nash ] = (/i + a)(l - (Si)) - 
(a + ji + 7 )(/i) < (// + a)(l — (Si)). Assume Ii(t) 7 ^ 0 for some t, and 
consequently I\ (t) 7 ^ 0 for all t. Therefore 

cm + T\ [pNash] (t) _ 7 + /X , /{[pNash](t) 

l|PN “ hK ' - -JW r ' 

and finally {Si[pNash]) = (7 + ^)(/V 1 ) = (So[p°pt]) and E[p opt ] = (fi + n)(l - 

(So)). □ 

Corollary 3. Let p^ as h and p opt he the Nash and optimal vaccination strate¬ 
gies associated to a given transmission parameter [3, respectively. Assume 
further that conditions in Theorem [3| and [3] are satisfied. Then E[p opt ] — 
E[p N ash] = r (7 + p + a)^- 1 ) > 0. 

Proof. From Equation 0, we have E[p opt ] = (p op tS 0 ) = (^ + a;)(l — (S 0 )), 
where (So) = (7 + p)(/3~). Furthermore, 

E[p N ash] = (PNash^l) = (jU + «)(1 - (Si)) - (/3-ZlSi) - a (If) . 

We also have that 

(Si) = ^ + ^ p ^ = ~Ui) + (7 +a 0(/2 x ) = (7 + v)(P 1 )- 

On the other hand, (j3SiIi) = r (7 + p)(/3^ 1 ). Furthermore, (Ji) = r(/3 _1 ). 
Finally, 

E[p Na sh] = (p + a) (1 - (7 + /i)(/3 -1 )) - r (7 + /^(/T 1 ) - ar(P~ l ) 

= E[p opt ] - r (7 + /i + «)(^ _1 ) ■ 


□ 
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In this work, instead of minimizing the vaccination effort in the set of 
strategies with TZq below one, we minimized the vaccination effort in a subset 
X P , which we have called the set of preventive strategies. Still, we can prove 
the following result. 

Proposition 2. TZo[p op t\ < 1. 

Proof. The inequality follows immediately from the Definitions [3] and [4] and 
Equation ([2]). □ 


4 Examples 


4.1 The constant case 

Let us consider /3 ft) = /3 0 , for all t. 

If /3o < 7 + /i, we have that Soft) = 1 for all t, and therefore {po := 
pft) = 0,Vf} € x P . As E[0] = 0, we conclude that Poptif) = 0. As there is no 
endemic solution, we conclude that PNash(t) = 0 . 

Now, assume (3 0 > 7 + p and assume additionally that r > 0 is small. 
From Theorem 2, p opt = (p + a)(S , 0 [p O pt]^ 1 - 1 ) = (// + cu) - lj which 

coincides with the optimal strategy in the traditional sense of lZo\p op t] = 1- 
The first condition on Theorem [3] is trivially satisfied and the second one is 
satisfied whenever 


r < r* := 


(p + «) (7 + P) ( Po 


7 + p + a \7 + p 


In this case, I\ [pN as b] = ^ 7 ^ 0, and therefore 0 = /i[pNash] / = 

A[PNash](/3o5'i[PNash] - 7 - p) and therefore Ai[p Nash ] = Furthermore 

0 S \ [pNash ft) p OL O' 1 1 [pNa-sh /^oh [PNash]‘S'l [PNash] P N as h S ] fp |\t as ] 

(/I + a)5i[PNash]- Finally, 


PNash (/xTOf) 


Po 


1 + P 


1 + 


a 


7 + P 


= Popt-r 1 + 


a 


7 + P 


"P Popt 


i.e., the rational level of vaccination will not be able to eliminate the disease 
as previously shown in (Bauch and Earn, 2004). It is clear that both p opt 
and PNash are admissible. 
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4.2 The sinusoidal case 

For the sinusoidal case, with j3(t) = /3o(l + £cos t) we can provide precise 
results. First, note that the condition Q in Theorem [2] is 

- . . - . . / /3n(l + £ COS t) 

—Possmt > — (/x + a)(3 o(l + £ cos t) (-1 

V 7 + A 4 



or, equivalently 


esin t < a )0o ^ _|_ £cos ^ 2 _ + a)(l + £cos t) . 


7 + A 4 


If 


£ < (/x + a) 


A) 


— 11— £(/x + a) 


2/3 0 


— 1 


,7 + A 4 / V7 + A 4 

condition Q will be satisfied for every t. It follows that last equation is true 

if 


£ <£ 0 \ = 


2 (/i + a) A) 

7 + A 4 


+ 1 — ji — a\ (/x + a) 


-i 


A) 


7 + A 4 


Now, assume f3(t) = j3 o(l + £cos t) with e < £ 0 . Note that 


(A)[popt]) = 


7 + A 4 


-2t r 


dt 


7 + A 4 


-2tt 


dt 


7 + A 4 


2vr 7 0 (3{t) 2nf3 0 ,/ 0 1 + ecos t /3 0 \/l - £ 2 

showing that the vaccination effort for the optimal solution 


E[p op t] = (a 4 + a)(l - ('S'o)) 


decreases with the oscillation amplitude. Furthermore, 


Popt (A) 

= (a 4 + a) 
— (fx + a) 
= (/i + a) 


A)(l + £ cos A) 


7 + A 4 


-1 - 


£^o sin t 


A) 


7 + A 4 

A) 

7 + A 4 


— 1J + £ 
1 | + £ 


yS 0 (l + £COSt) 

(A 4 + a )A) 


cos 


7 + A 4 

V'(A 4 + oQ 2 A) + (A 4 + 7) 2 

(a 4 + 7 ) 


t — sin t^ + (9 (£ 2 ) 

cos(t + <p) + O (£ 2 ) 


? 
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where p = arctg ■ It is clear that (£-) = 0 and therefore (p op t ) is the 

same as in the constant case. However, the vaccination effort (popt^olPopt]) = 
(/i + oi) (1 — (>5[Popt])) = (/i + oi) ^1 — 3 ^j== j is strictly smaller than in 
the constant case. Furthermore, to first order in £, the optimal vaccination 
strategy p opt lags behind the transmission rate by a phase shift of p. In 
particular, if the birth/mortality rate is very high (and therefore there is 
an extremely fast renewal of susceptible individuals), i.e., p —>■ oo, then 
<p —> arctg A. This means that the optimal vaccination time shift will depend 
on the average transmission rate. In the more realistic case /x —>■ 0, i.e, when 
the renewal is low, then ip —> arctg . This case is of much higher interest 
in practice as it models the case when the disease time scales (contact rate, 
recovery time, temporary immunity) are much smaller than the typical time 
scale of one generation of the population. Assuming « 1 (i.e., immunity 
lasts for approximately TZo ~ a~ l years), then p « |, i.e., in a seasonal 
epidemic, the peak of the vaccination rate, p, should be approximately 1.5 
months before the transmission peak. Note that this is not the peak of the 
instantaneous number of vaccinations, pS, as illustrated in Figure |3j 

Now, we use Theorem [3] to obtain the Nash strategy. Condition Q can 
be rewritten 

s sin t 

—£sin t + O (e 2 ) = -< 7 + 11 < /3 0 + e(/3 0 cos t — sin t) + O (e 2 ) 

v ' 1 + £ cos t v 7 

— Po + 1 cos(t + (p)+0 (e 2 ) , 

for a certain <p. This is true if, for example, £<^ + p<(5q~ E\J /3q + 1 , in 
particular if /3 0 > 7 + /i, and £ is small enough. 

Equation ([7]) is equivalent to the non-negativeness of PNash, given by Equa¬ 
tion (J 8 ]) . Let us assume that the relative risk of the vaccination is low, i.e., 
r = 0(e). After some extensive, but straightforward calculations, we con- 
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elude that 


PNash [£](*) 


=(p + a) 


A) 


7 + P 


— 1 — r 


a 


7 + P 


+ 1 


+ E- 


A 


7 + P 


a/( 7 + p ) 2 + 1 ( ^ + a COS (t + V’) “ 7T Sill (^ + 

V7 + P Po 


— (/i + <a) 

+ £ 


A) 


7 + P 
1 


a 


7 + P 


+ 1 


(7 + P) 2 


VTt + p ) 2 +1 \j Pi (p + «) 2 + (7 + p ) 2 cos(t + <p + VO, 


where tan tb = —j— and tan p = a ? + ^ - . It is clear that if -+ 2 - > 1, then, for 
e, r small enough, the Nash-equilibrium is given by PNash > 0 as obtained in 
the previous equation. It is also clear that both p opt and PNash are admissible 
at leading order. 

If /x —>• 0, A ) a /7 — 0(1), then tp « this means that pNash will peak 
shortly before p opt ; however amplitude oscillations are slightly larger for pNash- 
Finally, the difference between both vaccination efforts are given by 


E[p opt ] - E[p Nash ] 


r( 7 + p + a) 
27 t/3 0 


f 27T 


d£ 


/ 0 1 + ecost 


r( 7 + p, + a) 

A,vT^ 


Figure [3] shows the peaks and oscillations of the optimal and Nash strategies, 
for the sinusoidal case, and how the difference between them depends on the 
vaccination risk. 

In Figure [4j we consider three different vaccination profiles for sinusoidal 
transmission: optimal vaccination, p(t) = p 0 pt[P]] optimal vaccination in the 
case of average transmission rate, p = p op t[(/3)]; and optimal vaccination in 
the case of maximum transmission, p = p op t[P\ (used in Lemma [ 5 J. This 
example also illustrates the result of Corollary [lj namely, it shows that, 
for a given sinusoidal transmission rate /3, the effort associated to p 0 pt[{/3)} 
is larger than the effort associated to p op t [/3]- However, p opt [(/?)] does not 
prevent initial outbreaks. Additionally, we show that p op t [/3] prevents the 
initial outbreak but with a higher associated effort. 
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p(t) p(t)S(t.) 




Figure 3: Optimal and Nash strategies for the sinusoidal transmission with 
low and high vaccination risk, r. We use // = (80T)” 1 , 7 = 52/T, a = 2/T, 
with T = 2tt and consider j3(t) = 26(1 + 0.3cos(t)), implying = 1 and 


IZo = — = 7 T « 3. Left: Both p opt (blue, continuous) and pNash (red, dashed) 
oscillate in a synchronous way. The peak p opt is 7 t/4 before the peak of the 
transmission rate; the peak of pNash is slightly before. Right: Time dependent 
vaccination effort, in the two cases. Above: small vaccination risk, r = 0.005; 
PNash is higher than p opt in the beginning of the epidemic season and lower 
otherwise. Below: high vaccination risk, r = 0.1; p opt > PNash for all t, 
but the difference is larger when f3 is decreasing. For simplicity, we plot in 
all cases j3(t ) in gray dotted line (out of scale). The choice of parameters 
is inspired by influenza like epidemiological parameters: life expectancy of 
80 years, recovery time of 1 week, and temporary immunity of 6 months, 
and TZ 0 ~ 3 (cf. (Goeyvaerts et al. 2015)). I 11 the absence of vaccination, 
there is only one stable attractor of the solution ( S(t),I(t )) with period T 
(cf. (Kuznetsov and Piccardi, 1994)). 
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Figure 4: Prevention of initial outbreaks for different vaccination profiles. 
We consider all parameters, including the transmission rate, as in Figure [3] 
and assume initial conditions given by <So[p](0) —0.01 and 1(0) = 0.01, where 
So[p] (0) is given by Lemma[4j We consider three different vaccination prohles: 
P = Popt[P} (blue, continuous), p = p op t[{/3)] (red, dashed) and p = p op t[P] 
(green, dash-dot). Note that the optimal strategy in the case of the average 
transmission rate is unable to prevent initial outbreaks (however, /(f) —y 0 
when t —y oo). Vaccination efforts E[p opt [5]] = 0.1165 > E[p opt [(/3)] = 
0.1098 > E[p opt [/3]] = 0.1090, respectively. 
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4.3 A critical case 


Condition Q in Theorem [2j provides a lower bound on the derivative of 
the transmission coefficient f3. In particular, if [3 is increasing, Theorem [2] 
provides (at least, in principle) one optimal strategy. However, to apply 
Theorem |2j it is important that j3 does not decrease instantaneously. In this 
section, we will study a critical example for optimal vaccinations, where f3 
satisfies the critical condition f3' = — (p + a)(3 — lj in (0, 27 t) and is not 

differentiable at t — 0. As usual, f3 is periodic in M. Explicitly, we consider 


m 


_7 + At_ 

1 — K exp(—(/i + a)t)) ’ 


t e [0, 27t) , 


(9) 


where K e (0,1) is a constant. 

We consider a mollification of (3, /3 e , such that f3 e is differentiable, satisfies 
the condition Q, and f3 e —>■ (3 when e — > 0 pointwise. Let p e := p op t[/3 e \. It 
is clear that p £ {t) —> 0 for t ^ 0, and 


r -2 it 


p e dt —* (/i + a)2n 
= 27t(/x + a) 
= log 


7 + ^ 
1 


{P)~ 1 


27t(/x + a 

e 2(/i+a)i 


e 2(fJ,+a)ir _ 

log-1 

& 1 - K 


= log 


1 - K 

l — J{ e ~2(/J.+a)n 

1 - K 


p 2(n+a)n _ jy- 

2n (h + a) = log---—-h loge _27r(M+Q ° 


1- K 


= : T. 


Thus p e —* T Pop t- From Equation (jsj) , it is clear that for t G 

(0, 27 r) , p opt [lim f3 e ] (t) = limpopt [/3 e ](f), arid therefore p op t[lLn [3 e ] = lim p opt [f3 £ ] 
as measures. This shows that discontinuities in (3 will be associated to peak 
vaccinations. From Lemma [4j we have that 


%opt](0) 


(/i + a) / 0 T e (^ +q L t 2 

X — g— r—( p+o)t 


fe(-7d s 


e r / 2 (l — e h t+c 9 T ) 

1 — g-r ~{n+a)T ’ 


where 9(s) = 0 for s < 0 and 9(s) = 1 for s > 0 is the Heaviside function. 
Furthermore 


1 _ e -(v+°9 T 

S\Pe](t)p £ (t)dt -> Te- /2 ~ _ c _ r , (At+a)r < (/i + a)T , 
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Figure 5: Optimal strategy for the critical case, with relaxed (3{t) Left: Re¬ 
laxed version of [3{t ) given by Equation (|9| with K = 0.862, and the other 
parameters as in Figure [lj with T = 2tt. We consider a sequence f3 e of differ¬ 
entiable functions, such that f3 £ {t) = j3(t ) for t G (e, 27t — e) and the points 
27 t — £ and 2n + e = e are connected by a third order polynomial. Right: a 
sequence of p £ := p Q pt [/3 £ ]; note that ( Pe) E>0 resembles a delta-sequence. 


where we used that 


xe 


-x/2 


l-e~ y 


= x 


sinli v, 

< y , Vx, y > 0 . 


sinh 


1 - e~ x ~y 

This last inequality follows from the convexity of the function 


sinh a = sinh 


^ ;0 H- : —-(a +/U)''] < —sinh0-1 -— sinh (a. + /3) , 


ol T /3 ol T /3 


ol T f5 


ol T f3 


and therefore, for a, (3 > 0, 

sinh a a 

sinh(a + (3) ~ a + (3 ~ f3 

Optimal strategy for the critical case, using a relaxed version of the trans¬ 
mission function {3(t ), is illustrated in Figure |5] 


5 Discussion 

In this work, we consider a SIRS model with periodic transmission, where 
we introduce periodic vaccination of adults. We are naturally led to consider 
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temporary immunity, as discussed in (Onyango and Muller 2014): in the 


case of diseases with life-long immunity and lifespan much longer than the 
period (// ~ a <C 1/T) the vaccination is not affected by periodicity. For 
instance, measles has periodic transmission rate but, since it confers perma¬ 
nent immunity, periodic vaccination is not used. In this work we implicitly 
assume // <C 1/T ~ a. 

We study the consequences of two extreme types of vaccination strategies: 
mandatory vaccination, where a certain predefined fraction of the population 
is vaccinated; and voluntary vaccination, where individuals can choose freely 
to be vaccinated or not, according to their risk perception. Classically, the 
objective is to minimize the vaccination effort while reducing the effective 
reproductive number below one, which guaranties long term disease elimi¬ 
nation. Here, we choose to work with an alternative definition of optimal 
vaccination. We define a class of preventive vaccination strategies as vac¬ 
cination profiles that, for any sufficiently small perturbation of the disease 
free state, the number of infectious individuals is monotonically decreasing, 
avoiding the occurrence of any epidemic event. This approach allows, for 
specific regular transmission functions 0(t), the derivation of an analytical 
expression of the optimal strategy. In general, we prove the existence of an 
optimal strategy, in a suitably defined closure of the space of all preventive 
strategies, which minimizes the vaccination effort. 

In this work, we extend the classical results by Bauch and Earn (2004) 
to periodic functions, based on a series of recent results on periodic diseases. 
We model human behaviour using classical economical theory, where indi¬ 
viduals are assumed to be rational and fully informed. We define the set 
of vaccination strategies that provide herd immunity, for which the rational 
strategy of a given focal individual is not to be vaccinated. Finally, we prove 
the existence of a Nash vaccination strategy as the strategy that minimizes 
the joint risk for every individual, taking into account the strategy of all 
other individuals. 

In general, both optimal and Nash strategies will not be functions but 
Radon measures. For specific forms of the transmission rate, we provide 
explicit formulas, which includes some important examples as constant or 
sinusoidal transmission functions. 

There are several natural limitations of the work presented here. One first 
limitation is that we consider only the stationary solution of the System (Jl|, 
but we never discuss the approach to this equilibrium. This is an important 
question, both in the study of ordinary differential equations (i.e., the study 
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of the basin of attraction) and in evolutionary game theory, where the study 
of Co'-limits of conveniently defined dynamical equations is preferred to the 
static study of Nash-equilibria. In the non-stationary case, a rational decision 
will require the ability to forecast the evolution of the epidemic, i.e., rational 
decisions will depend on future decisions of the entire population and not 
only on the past decisions. This is mathematically described by the so called 
“mean field game theory” (Lasry and Lions, 2007) and will be object of a 
future work. 

Closely related ideas will also help us to solve one of the major gaps of 
the current work: the lack of a numerical method for finding Nash-equilibria 
solution when Theorem [3] fails. More precisely, the idea will be to develop a 
numerical method that allows constant update in individual decisions and, 
consequently, also at the population level. As discussed before, this will 
require, at the individual level, a certain expectation on the future evolution 
of the disease. 

It is also important, and will be subject of a future work, to design a 
precise scheme, possibly numerical, that allows to go beyond Theorem [2j 
This will require the use of Optimal Control Theory. In fact, given p, it is 
possible to explicitly obtain the disease free solution S[p](t) (see Lemma [4| 
and therefore we need to minimize p(t)S[p](t)dt in y p [/3]. (Equivalently, 
we may maximize S\p](t)dt in the same set, as (pS 0 ) = (p+a) (1 — ((So)).) 
We also plan to compare p opt with the optimal solution in the approach in 
which the vaccination effort is minimized in the class of vaccination functions 
p such that TZq < L 

Furthermore, despite the simplicity of the periodic SIRS system (even 
with vaccination), solutions can be extremely complicated; even chaotic so¬ 


lutions may be present in such simple systems, cf. (Kuznetsov and Piccardi 


1994). The coupling of the differential equations with human rational be¬ 


haviour presented in this work only started the exploration of all this math¬ 
ematical richness. 


A Proof of Lemma ED 


We follow closely the proof at (Rebelo et ah, 2012), where x\ — /, x% — S 
and x% = R. Also, m — 1 indicates that there is only one infectious class 
and n = 3 denotes the three possible classes in the model. We readily verify 


that conditions (Ai) — ( A 5 ) in (Rebelo et ah, 2012) are satisfied. Unique 
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ness and stability (in the disease free subspace) of the disease free solution 
(So(t), 0, Ro(t)) is guaranteed by standard theorems. The linearisation of sys¬ 
tem (|TJ) restricted to / = 0 around the disease-free solution (So(t), 0, Ro(t)) 
is given by 

j s' = —fis — p(t)s + ar 
1 r' — —fir + p(t)s — ar 

that can be explicitly solved to get 

( s(t ) = s(0)a(f) + (s(0) + r(0))6(t) 

\ r(t) = (s(0) + r(0))e _/it — s(t) 

where a(t) = e ~ S Q r+a+pir) dr anc [ _ j* ae -/*T-f*(ti+a+p(i))di This 
yields that the monodromy matrix M(t) of the linearised system is 


T a(T) + b(T) b(T) 

[ ' [ e~^ T - a(T ) - b(T) e~» T - b{T ) ' 

We compute the Floquet multipliers p\ = a(T ) < 1 and p 2 = e _AlT < 1 and 
conclude condition (Ag). We verify immediately that conditions (A?) and 
(A 8 ) are also satished. 

Let (S,I,R) be a solution of the system (Jl| and (So, 0, R 0 ) the disease 
free solution of the same system. Therefore 


(S-S 0 y = - aI-/3IS-(p(t)+fi+a)(S-S 0 ) < -(p(t)+fi+a)(S-S 0 ). (10) 


This yields for every t > 0. Consequently, by Gronwall’s lemma 

S(t) - S 0 (t) < (. S(t 0 ) - So(i 0 )) exp £ ( p(s) + fi + a)ds 

for any t > t 0 > 0. 

So, for any e > 0 there is ^(e) > 0 such that for any t > ti we have 

S(t ) - S 0 (t) < e. (11) 

Now, assume that there is to > 0 such that I(t) < e for every t > to. 
Therefore, as S < 1, 

(So-S)' = aI+piS-(p(t)+ii+a)(S Q -S) < (0+a)e-(p(t)+fi+a)(S Q -S), 

(12) 
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and, by Rebelo et al. (2012, lemma 1) there exists k > 0, independent of e, 
and t 2 (e) > t 0 such that for all t >t 2 


S 0 (t ) - S(t) < ke. 


(13) 


For any solution in the disease-free subspace (i.e., with I(t) = 0 for all 
t > 0), we have the validity of Conditions ( fToj ) and ( Jl2j) , and, therefore, we 
conclude that the disease free solution is globally asymptotically stable. 

Now, we show the alternative in Lemma 0 Let IZn = ' :l3So) , as defined 


m 


Onyango and Muller (2014) and in a more general setting in Rebelo et al 


(2012); Wang and Zhao (2008). We show that the only relevant assumption 
in |Rebelo et al. (2012 theorem 2) is the value of TZq. In particular, we define 
the lxl matrices F(t ) = [f3(t)So(t)] and V(t) = [/i + 7 ] (Rebelo et al. 2012). 


By (11), for every e > 0 there is t x > 0 such that for t > 1 i, 


i\t) = mms(t) - ( 7 + m o < + s 0 m - ( 7 + 

* (W> - m ) m 

where Ai(e) = mirp ■ Notice that Ai(e) > 0 and that Ai(e) —* 1 from 
below when e —> 0. 


If IZq < 1, Rebelo et al. (2012, theorem 2, condition 1) guarantees that 


the diseases dies out, I(t) —> 0 as t —> oo and that (S'o,0,i?o) is globally 
asymptotically stable. 

Now, assume TZ 0 > 1; from A /3(t)S 0 (t) dt — (n + 7 ) > 0 and by 
continuity in t, we have that F(t) — V(t) is irreducible for some t G [0, T ). 
For any e > 0, if there is t 0 > 0 such that I(t) < e for every t > t 0 then, 


by (13), there is t 2 > to such that, for t > t 2 , 


m > mm(Mt) - m - p + -dm > m, 

where A 2 (e) = max t So(t)-ke sa fisfi es li m e-s>o+ A 2 (e) = 1 and A 2 : (0,e*) —>■ M + 
if we choose e* < ^ mirp So(t) (observe that Lemma[4]and its proof guarantees 
that min f S'o(t) > 0). 

We conclude that the conditions in Rebelo et al. (2012| theorem 2, state¬ 
ment 2) are satisfied and there is uniform persistence of system (Jl| with 
respect to /. 
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Finally, we prove the existence of a persistent periodic solution. We define 
the T-mapping P : A 2 —» A 2 by P(x o) = x(T, (cco, 0)), where x(-, (xo,0)) is 
the solution of System [l] with initial conditions x(0) = Xo E A 2 := {( S, I ) G 
S + / < 1}, the two-dimensional simplex. P is a continuous map such 
that P(M 0 ) C M 0 for M 0 = {(I,S) E A 2 : / 7 ^ 0}. Observe that M 0 is 
an open set of A 2 with the topology induced in A 2 . As we have uniform 
persistence of system (|TJ) with respect to / we also have uniform persistence 
of P with respect to Mo as described in Zhao (1995) (for a more general 
case see ( Magal and Zhao] , 2005)). Applying (Zhao, 1995, theorem 2.1) we 
conclude that P : M 0 —> M 0 admits a global attractor and there is a fixed 
point for P in that attractor, which is a T-periodic solution of (|T|) (see, for 
example, lemma 4.4 in Verhulst (1996)). 


B Proof of Theorem U 


Proof. First we recall that RM + is compact with the weak topology, cf. (Ko- 
ralov and Sinai 2007). We have that C([0,T]) C RM + in the sense that 
for each continuous function we consider the correspondent cumulative dis¬ 
tribution function. Consequently the closure C([0, T]) C RM + is compact. 
The map E : C([0,T]) —> [0,+ 00 ] defined in Definition [l] and extended in 
Lemma [3] is continuous with respect to the weak topology. 

Existence of p Q pt : As xf C RM + is compact, from the continuity of E, 
we conclude that there is a measure p opt G xf suc h that for all p E y p , 
E[p] > E[p opt ]. 

Existence of pNash : Let us consider hxed time intervals At, such that 
T / (At) = N E N, and consider periodic continuous piecewise affine func¬ 
tions in intervals (iAt, (i + l)At), i E {0,..., N~ 1} (i.e, functions, such that 
/(f) = f([t/At\ At) + ■d(h/ A d+ 1 )^)-/(h/ A d A <) (j. _ (t/At\At)] furthermore, 
f(NAt) = /(0)). These functions can be represented by vectors in Let 
v G T := {v G M^|0 < v± < v 2 < • • • < Vn < (a + p)NA^}. The set T is 
convex and compact. Now for each vector v G T consider the distributions 
p v : [0, T] —>■ M + , such that the associated cummulative distribution is given 


by P v (f) = Vi + 


v i + l-Vi 

At 


(t — iAt) for t E [iAt, (i + l)Af], i E {0,..., N — 1}. 


Consider I[p v ] and S\p v ] solutions of System (JTJ) , and define for v,v* G T 
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the joint risk (except for some immaterial constants) 

N 

P A V,v] = -^{(r - (3(iAt)I[p v ](iAt)) e~ v >} . 

i =1 

It is clear that v* « J 0 * At dP* = P*([ 0 , iAt]). We define a function T : 
T —* 2 t such that v G J-’fv] if and only if p A< [v, v] < p At [v', v] for all v' G T. 

It is clear that W[v] ^ 0, as T is compact and p At [-,v] is continuous. 
Now, we prove that J-"[v] is closed and convex. The first property follows 
again from the continuity of p At [-, v]. Define r* = r — /3(iAt)I[p v \(iAt). We 
divide the last property in two cases: 

1. Assume that A < 0 for all i and let v* be such that p[v*,v] < p[v',v] 
for all v' G T. Assume in addition that there is v such that p[v*, v] = 
p[v, v] for all v G T. Therefore, for a G (0,1), 

p At [aw* + (1 - a)v, v] = - Ae-“<- (1 - a)ii 

l 

< — fj ( ae~ Vi + (1 — a)e _i;i ) 

i 

= ap At [v*, v] + (1 - «)p Ai [v, v] = p[v*,v] 

< p At [aw* + (1 - cc)v, v] ; 

hence v*, v G J-’fv] implies that av* + (1 — a)v G T . 

2. Let X := {z|f; > 0} 7 ^ 0. In order to minimize p At [-,v], we impose to 
each j G I the minimum possible value, i.e., v* = v*_ x . Therefore, we 
shall minimize 

p^[v*,v] := -^{(r-/3(iAt)/[p v ](iAf))e-<} . 

The existence of a minimum is guaranteed by the compacity of Tx := 
{v G M A |0 < v 1 < v 2 < ■ ■ ■ < v N < (a + p)N^, i e Z => ry_i = 1 ^}. 
Then, we repeat the previous analysis and conclude that X”[v] is closed 
and convex. 

We conclude that the set of best replies is non-empty, convex, closed 
and due to the continuity of I in p (see Lemma [lj and of p At in v and v*, 


34 


the graph is closed. Therefore from standard applications of Kakutani fixed 
point theorem, there is a fixed point vector of the function tF, v (A ^, such 
that its affine function continuation p^ A d is a Nash-equilibrium restricted to 
affine functions with steps At. See, e.g., 

Furthermore 

7V-i TV-1 5 a 

p( At )([0, T]) = p( At ^(zA£)At < ViAt < (a+p)N ——— At = (a+p)T———- 


From the compactness of C([0, T]), there is a measure p G C([0, T]) such that 
liniAt-KjP^ = P (possibly after taking subsequences), where the convergence 
is in the weak topology. 

The last step is to prove that p is indeed a Nash-equilibrium in C([0, T]). 
Assume it is not; then, there is p in C([0, T]) such that p\p,p] < p[p,p]. 
From the continuity of p, there is a (A£)o > 0, small enough, such that all 
restricted Nash-equilibria found above are such that p[p, p^ A ^] < p[p^ At \ p^] 
for At < (At) 0 . Let p n be a sequence of continuous functions in (7([0,T]) 
such that p n —» p weakly, and therefore p[p n ,p^ At ^] < p[p( A b,p( A d], f or a 
certain value of At and n large enough. Using the fact that p At [v*, v] is the 
trapezoidal approximation of p\p*,p\ (and therefore differs in 0((At) 2 )) and 
taking At possibly even smaller, we conclude that p( A d i§ not a restricted 
Nash-equilibrium, contradiction. □ 


(Osborne and Rubinstein, 1995). 
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